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We investigate the Cosmic Censorship Conjecture by means of the horizon wave-function (HWF) 
formalism. We consider a charged massive particle whose quantum mechanical state is represented 
by a spherically symmetric Gaussian wave-function, and restrict our attention to the superxtremal 
case (with charge-to-mass ratio a > 1), which is the prototype of a naked singularity in the classical 
theory. We find that one can still obtain a normalisable HWF for a 2 < 2, and this configuration has 
a non-vanishing probability of being a black hole, thus extending the classically allowed region for 
a charged black hole. However, the HWF is not normalisable for a 2 > 2, and the uncertainty in the 
location of the horizon blows up at a 2 = 2, signalling that such an object is no more well-defined. 

This perhaps implies that a quantum Cosmic Censorhip might be conjectured by stating that no 
black holes with charge-to-mass ratio greater than a critical value (of the order of y2) can exist. 

PACS numbers: 04.70.Dy,04.70.-s,04.60.-m 


I. INTRODUCTION 

A complete understanding of the gravitational collapse 
of a compact object remains one of the most challeng¬ 
ing issues in contemporary theoretical physics. The gen¬ 
eral relativistic (GR) description, resulting in the forma¬ 
tion of a black hole (BH) or naked singularity (NS), was 
first investigated in the papers of Oppenheimer and co¬ 
workers [1]. Although the literature on the subject has 
grown immensely (see, e.g. Ref. [2]), many technical and 
conceptual issues remain. One of these is the famous Cos¬ 
mic Censorship Conjecture (CCC), proposed by Penrose 
in 1969 [3], which states that no singularities will ever 
become visible to an outer observer in a generic grav¬ 
itational collapse starting from reasonable nonsingular 
initial states. To date, the conjecture remains unproved, 
and it is considered one of the most important open prob¬ 
lems in gravitational physics. Another great open issue in 
GR is the problem of considering the quantum mechani¬ 
cal (QM) nature of the collapsing matter [4]. We will here 
address both issues for the Reissner-Nordstrom (RN) ge¬ 
ometry, which describes charged BHs, a subject of many 
theoretical investigations in the past (see, e.g. Ref. [5]). 

Most attempts at quantising BH metrics consider the 
gravitational degrees of freedom unrelated to the matter 
state that sources the geometry. More recently, the Hori¬ 
zon Wave Function (HWF) formalism was proposed [6], 
as a way of quantising the Einstein equation that deter¬ 
mines the gravitational radius of a spherically symmetric 
matter source and its time evolution [7], which instead 
relates the quantum state of the horizon to the quan¬ 
tum state of matter. This formalism was then applied 
to a few different case studies [8-10], yielding apparently 


sensible results in agreement with (semi)classical expec¬ 
tations, and there is therefore hope that it will facilitate 
our understanding of the formation of BHs from QM par¬ 
ticles. In particular, it seems natural to extend this for¬ 
malism beyond Schwarzschild BHs and tackle the CCC 
from a quantum perspective by considering an electrically 
charged particle represented by a Gaussian wave-packet 
in the classical regime in which it would be a NS. 


II. ELECTRICALLY CHARGED SPHERICAL 
SOURCES 


We start by recalling the classical RN metric can be 
written as 


ds 2 

with 


—/ dt 2 + / _1 dr 2 + r 2 (d0 2 + sin 2 0 defy 2 ) , (1) 


/ 


2£ P M Q 2 

"F 2 5 

m p r 


( 2 ) 


where M and Q respectively represent the ADM mass 
and charge of the source, £ p is the Planck length and m p 
the Planck mass 1 . For \Q\ < £ p M/rn p , the above metric 
contains two horizons, namely 


R± = 4 — ± 



Q 2 , 


( 3 ) 


and represents a BH. The two horizons overlap for 
| Q | = £pM/m p , the so-called extremal BH case, while 
for |Q | > £ p M/rrip no horizon exists and the central sin¬ 
gularity is therefore accessible to outer observers. This 
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1 We shall use units with c = kg = 1, and always display the 
Newton constant G = £ p /m p , so that h = i p m p . 
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is the prototype of a NS, which we will refer to as the 
“superextremal geometry”. It is in fact more convenient 
to express all relevant quantities in terms of the mass M 
and the (positive definite) specific charge 


M 

£ p M 


m p 


( 4 ) 


Using this parameter, the above expression (3) becomes 


R±=£ p — (l±y/l- a 2 ) , (5) 

and the three regimes mentioned above are then explic¬ 
itly parametrised as i) 0 < a < 1 for the BH with two 
horizons , ii) a = 1 for the extremal BH, and iii) a > 1 
for the superextremal geometry. 

We shall now investigate the superextremal geometry 
from a quantum mechanical perspective by first deter¬ 
mining the HWF for a < 1 and then extending it contin¬ 
uously into the regime a > 1. 


A. HWF for Gaussian sources 

The general procedure that leads to the HWF [6, 7] is 
based on lifting the gravitational radius Rh of a spher¬ 
ically symmetric QM system to the rank of a quantum 
operator. This step can be physically motivated by first 
recalling that the coordinate r in a spherical metric like 
the one in Eq. (1) is invariantly related to the geometrical 
area 4 it r 2 of a sphere centred on the origin r = 0, and 
is therefore a natural candidate to become an observable 
in the quantum theory. Moreover, the specific property 
that qualifies r = R]j is that it represents the location of 
trapping surfaces and thus determines the causal struc¬ 
ture of the space-time, which one can assume will also 
ought to remain an observable property in the quantum 
theory. In details, we recall that in a neutral spherically 
symmetric system, Rn(r) = 2£ p M(r)/m p , where 

M(r, t) = ( p(r,t)f 2 dr , (6) 

Jo 

is the Misner-Sharp mass. One should notice that M 
represents the total energy (thus, roughly speaking, in¬ 
cluding the negative gravitational energy) and is related 
to the energy density p of the source via the flat space 
volume. A specific value of r is then a trapping surface if 
Ru(r) = r, whereas, if Rn(r) < r, the gravitational ra¬ 
dius is still well-defined but does not correspond to any 
peculiar causal surface. In the electrically charged case, 
we have two gravitational radii and corresponding oper¬ 
ators, namely R±. The classical relation (5) will then be 


2 We remark in passing that the inner horizon gives rise to an 
instability usually referred to as “mass inflation”, but we have 
investigated under which conditions the inner horizon is actually 
realised in the quantum context in a separate work [11]. 


reinterpreted in this context as the operatorial equation 
relating R± to the total energy M of the system, with 
R± acting multiplicatively on the HWF and M acting 
multiplicatively on energy eigenstates. Finally, we will 
consider the ratio a as a simple parameter. 

Let us consider as a source of the RN space-time an 
electrically charged massive particle at rest in the ori¬ 
gin of the reference frame, represented by a spherically 
symmetric Gaussian wave-function 

r 2 

~~ £ 3/2 ^3/4 ' 

We emphasise that the radial variable r in the above is 
interpreted to be the same as the coordinate r in the 
metric (1), and is therefore a measure of the particle’s 
size as seen from an outer observer who, e.g. scatters 
particles against it. We shall also assume that the width 
of the Gaussian t is the minimum compatible with the 
Heisenberg uncertainty principle, that is 3 

777 

i = A m ~ l p -*■ , (8) 

m 

where A m is the Compton length of the particle of rest 
mass in. In momentum space, the wave-function of the 
particle described above is 

_ p 2 ^ 

^(P) = A e 3/C 3 /4 ’ ( 9 ) 

where p 2 = p ■ p is the square modulus of the spatial 
momentum, and the width A = m p £ p /£ — m. For the 
energy of the particle, we shall employ the relativistic 
mass-shell relation in flat space, 

M 2 = p 2 + m 2 , ( 10 ) 


in analogy with the expression of the Misner-Sharp 
mass (6). This choice does not therefore imply that the 
effects of curved space-time are a priori discarded, but 
rather that they should be included in the very definition 
of total energy of the system 4 . 

For a < 1, it is clear that one can now write a HWF 
for each of the two horizons. In fact, from the quantum 
version of Eq. (3), the total energy M can be expressed 
in terms of the two horizon radii as 


and 



R+ +R- 
2 


( 11 ) 


R± 


— 


1 ± Vl - a 2 

1 =F \/l - a 2 


( 12 ) 


3 A thorough analysis of the changes that occur by relaxing this 
condition can be found in Ref. [7]. 

4 Let us also remark that, for a source of Planckian mass and size, 
quantum effects might make the very concept of “curved space- 
time” inadequate in its proximity, whereas determining (asymp¬ 
totic) cross sections for particle collisions might still be meaning¬ 
ful. 
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Note that we promoted M, R+, and P_ into operators 
M, R+, and P_, which are related to the corresponding 
observables. Our specific choice is not unique, and it 
is associated with usual ambiguities when going from a 
classical to quantum formalism. 

The unnornralised HWFs for R + and R- are then ob¬ 
tained by expressing p from the mass-shell relation (10) 
in terms of the energy M in Eq. (11), and then replacing 
one of the relations in Eq. (12) into Eq. (9). The two 
HWFs are then given by 


lfal(R±) = A/± 0 (R± - -Rmini) 


x exp 


Wp R 


2 

± 


2 A 2 £2 (1 ± \/T 



,(13) 


where the Heaviside function arises from the mini¬ 
mum energy in the spectral decomposition of the wave- 
function (7) being M = m , which corresponds to 


Rm m + = £p — (l ± \/l — a 2 ) • (14) 


Finally, the normalizations A f± are fixed by assuming the 
scalar product 5 


W’h I M = 4tt / ^(R±)MR±)R±dR± , (is) 
Jo 

where, like in the previous equation, the upper signs are 
used for the normalization of ipn (P + ), while the lower 
signs are used when normalizing f/>H(P-)- 

The probability density that the particle lies inside its 
horizon of radius r = R± can now be calculated starting 
from the wave-functions (13) associated with (7) as 


V<±(r<R±)=Ps(r<R ± )Vn(R ± ) , (16) 

where 

rR± 

Ps(r < R±) = 47r / \i/>s(r)\ 2 r 2 dr (17) 

Jo 

is the probability that the particle is inside a sphere of 
radius r = P±, and 

Vh(R±)=4ttR 2 ± \^(R ± )\ 2 (18) 


is the probability density that the sphere of radius r = 
R± is a horizon. Finally, one can integrate (16) over all 
possible values of the horizon radius R+ to find the prob¬ 
ability for the particle described by the wave-function (7) 
to be a BH, namely 

poo 

Pbh+ = / V <+ (r < R+) d R+ . (19) 

J ^min + 


5 The analytical expressions for the normalization of the HWFs 
are very cumbersome and not particularly significant, thus we 
will omit them throughout the paper. 


The analogous quantity for P_, 

poo 

Pbh- = / P<-(r < R-)<\R- , (20) 

'' Rmin — 

will instead be the probability that the particle lies fur¬ 
ther inside its inner horizon, and both and R + are 
therefore realized (for more details about this case, see 
Ref. [11]). 

Before moving to the specific topic of this work, let 
us note that for a particle with mass M — m smaller 
than the Planck scale, the classical values of R± from 
Eq. (5) are smaller than the Planck length. This how¬ 
ever should not prevent one from considering the cor¬ 
responding R± as valid quantum operators representing 
the gravitational radii of the system. In fact, let us first 
point out that, for any value of the particle’s mass, if 
( R ±) < £, these radii most likely do not correspond to 
horizons, exactly like in the classical theory there is no 
trapping surface of size r if the corresponding gravita¬ 
tional radius R(r) = 2 £ p M(r)/m v < r. Consequently, 
the probability that the particle is a BH decreases very 
fast below the Planck mass, as it was shown for the sim¬ 
pler case of an electrically neutral particle in Refs. [6-8], 
and is confirmed for the probability Pbh+ of the RN BH 
in Ref. [11]. This means that, although one can have 
(R± ) < £ p , QM fluctuations due to the uncertainty in 
the particle’s size dominate below the Planck energy and 
the horizon is simply not realised as an actual trapping 
surface (see also footnote 4). More realistically, it was 
shown in Ref. [9], that a BH has a significant proba¬ 
bility to form from the collision of sub-Planckiian par¬ 
ticles only if the center-of-mass energy reaches into the 
trans-Planckian regime. Finally, in Ref. [8], it was shown 
that the uncertainty in the horizon size obtained from 
the HWF, combined with the standard Heisenberg uncer¬ 
tainty for the particle’s size, yields a minimum detectable 
length which is always larger than the Planck scale. In 
other words, a minimum detectable length of the order 
of the Planck scale naturally emerges within the HWF 
formalism, without the need to assume it, and without a 
modification of canonical commutation rules. 


B. Superextremal geometry 

We will now focus on studying overcharged sources, 
represented by the range of specific charge a > 1. It 
is well known that in the classical theory of gravity, the 
CCC a priori forbids the existence of NSs. In the case 
of the classical charged BHs, this would precisely corre¬ 
spond to a > 1, so it is interesting to investigate whether 
quantum physics leads to any modifications or predictions 
therein. Our guiding principle will be to assume that the 
quantum states in the regime a > 1 can be obtained by 
extending continuously the HWF from the case a < 1. Of 
course, this is by no means the only possible choice, but 
one can at least hope that it leads to consistent predic- 
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tions for a not too much larger than the classical limiting 
value of 1. 

The first problem that we encounter for a > 1 is that 
the operators R± are not Hermitian. We could stop right 
there and say that there are no observables which corre¬ 
spond to R± in the regime a > 1. However, we will 
make a simple (and perhaps non-unique) choice, which 
will allow us to proceed and probe the classically for¬ 
bidden region. We will take the real parts of the mul¬ 
tiplicative operators R ±, which are certainly Hermitian, 
to correspond to quantum observables. With this choice, 
the modulus squared of the two HWFs are given from 
Eq. (13), for R ± > R min ±, by 


|V'h(^±)| 2 = A4 2 exp 


to 2 R 


2 

± 


A 2 / 2 (1 ± Vl - a 2 ) 2 


,( 21 ) 


which, for a > 1, becomes one expression 


IV'H {R)\ 2 = A/" 2 exp 


2 - a 2 m 2 v K 2 | 

A 2£fj ’ 


( 22 ) 


where R now replaces both R + and f?_. This HWF is 
still normalizable in the scalar product (15) if R belongs 
to the real axis and the specific charge lies in the range 

1 < a 2 < 2 . (23) 


We could therefore infer that no normalizable quantum 
state with a 2 > 2 is allowed, or that there is an obstruc¬ 
tion that prevents the system from crossing a 2 = 2. This 
point will be further clarified after we have fully deter¬ 
mined the HWF. 

The fact that the HWF (22) is the same for R + and 
i?_ reflects the classical behaviour according to which 
the two real horizons merge at the critical value of a = 1, 
and then mathematically extend as one into the complex 
realm for a > 1. We then need to address what happens 
to the Heaviside function in Eq. (13) when we extend it 
into the superextremal regime. First of all, we note that 
although Eq. (14) becomes complex for a > 1, its real 
part is again the same for R + and , namely 


Rn lin — Ho 


(i±vT 



(24) 


We can then show that the same continuity principle, 
which led us to Eq. (22), requires that R be bounded 
from below by this R m in- In fact, the expectation value 
for R is, in this case, 


[R) =45 


IV’h (R)\ 2 R 3 dR 
2 £\ (2 — a 2 + a 4 ) /£ 


(25) 


a/2 — a 2 2 a/ 2 — a 2 + a 2 e <* 4 ypK erfc ^ ^ 2 a2 a2 ) 



FIG. 1. The expectation value ( R) (solid line) and its uncer¬ 
tainty A R (dashed line) as functions of the specific charge for 
1 < a 2 < 2 and m = m p (£ = £ p ). 


and it matches exactly the corresponding expressions for 

a < 1 [11], 


[ R± ) = 4 7r 


IV’h(H±)|“ R± d R± 


4/ 2 (l ± Vl-a 2 ) /(. 
2 + e y/n erk(l) 


namely 


U 2 Ji 


lim (R) = 
a\i 2 + ey/Ti erfc(l) a/'l 


= lim ( R± ' 


One can likewise show that the uncertainty 

AR 2 (£, a > 1 ) = (R 2 ) - (R) 2 , 
matches the corresponding uncertainties 

AR 2 ± (£,a<l) = (R 2 ± )-(R ± ) 2 , 


(26) 

(27) 

(28) 

(29) 


at the specific charge a = 1, but we omit the explicit 
expressions since they are rather cumbersome. We just 
note that, for a = 1, the width of the Gaussian i > 
( R ) for to < A/2+^/^erfc(l) TOp/2 ~ 0.8 to p , so that 
quantum fluctuations in the source’s size will dominate 
for masses significantly smaller than the Planck scale (in 
qualitative agreement with the neutral case [6, 7, 11]). 

It is now interesting to analyse the limit a 2 —> 2. One 
may have already noticed that 


(R) 


2 5 / 4 £ 2 /£ 

\J 7t (a/ 2-a) 


(30) 


so that the ratio (R)/£ blows up at a 2 = 2 for any values 
of the mass m = m p £ v /£. The same indeed occurs to the 
uncertainty, since 

AR ~ a/3tt/ 8-1(-R) ^ 0.4(H) , 
for a 2 —> 2 (see also Fig. 1). 


(31) 
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FIG. 2. Pbh as a function of a for m = 2m p (solid line), 
m = m p (dotted line) and to = 0.5 m p (dashed line). Cases 
with to S> m p are not plotted since they behave the same as 
to = 2 to p , i.e. an object with 1 < a 2 < 2 must be a BH. 


Using Eq. (19), one can also calculate the probability 
Pbh that the particle is a BH for a in the allowed su¬ 
perextremal range (23). This probability is displayed in 
Fig. 2. One notices that, for a particle mass above the 
Planck scale, Pbh is practically one throughout the en¬ 
tire range of a (thus extending a similar result that holds 
for a < 1 [11]). Moreover, even for to significantly less 
than to p , Pbh approaches one in the limit a 2 —> 2. We 
recall here that Pbh <C 1 for small to is essentially re¬ 
lated to £ (R), and the system is thus dominated by 

quantum fluctuations in the source’s position well below 
the Planck scale. On the other end, since both ( R ) and 
A R blow up on approaching a 2 = 2, the superextremal 
configurations with a significant probability of being BHs 
contain strong quantum fluctuations in the horizon’s size. 

Let us conclude this section by emphasising that, in or¬ 
der to achieve the above results, we had to choose a way 
to continue the HWF and operators that are straightfor¬ 
wardly defined for a < 1 into the overcharged regime. 
Since this choice is not apparently unique, one might 
wonder whether different options would lead to signifi¬ 
cantly different outcomes. Although we have not per¬ 
formed a full survey, it is important to remark that we 
were able to ensure the continuity of expectation values 
across a = 1, like in Eq. (27), by imposing simultane¬ 
ously that R equals the real part of R± and Eqs. (24) 
and (22), whereas we found no ways to make other smil¬ 
ingly natural options, like R = |P±|, work. 


can still make a quantum BH. The basic reason for this 
is that in our formalism the location of the horizon is not 
given by a sharp classical value, instead it is described by 
a quantum wave function with associated uncertainties. 
Second, the clrarge-to-mass ratio a cannot be arbitrarily 
large, even in the context of QM. We found that for a 2 > 
2 the HWF cannot be normalised, and thus it is not 
describing a well defined physical object. Moreover, at 
the same value of a 2 =2, the uncertainties in the location 
of the horizon become infinite, signalling again that such 
an object stops being well defined. We should warn the 
reader that the specific value of the upper limit a 2 = 2 
in Eq. (23) might simply be a consequence of describing 
the source as the Gaussian function (7), and should not 
be taken literally. However, it is likely that the overall 
qualitative picture remains in a more general context, 
and our results imply that perhaps a quantum version 
of the CCC might be formulated by stating that no BHs 
with the charge-to-mass ratio greater than a critical value 
(of order \/2) can exist. 

We should here recall that for the charged R.eissner- 
Nordstrom metric with a < 1 analysed in Ref. [11], as 
well for neutral sources [7, 8, 10], the single Gaussian con¬ 
stituent (7) leads to unacceptably large uncertainties in 
the horizon size of large astrophysical BHs. In fact, one 
has A R ~ (R), even for very large mass in, for which 
one expects a semiclassical behaviour for the horizon size 
R. The above quantum CCC will therefore have to be 
tested further, by considering models of BHs that allow 
for a semiclassical limit A R <C (R)- An example of 
such models is given by those in Refs. [10, 12], which 
contain a very large number N of light constituents, 
whose wave-functions span the entire region inside R, 
and A R/(R) ~ A -1 . The emerging picture is that BHs 
of any size should be treated as macroscopic quantum 
objects (just like superconductivity and superfluidity are 
macroscopic quantum phenomena at scales where one ex¬ 
pects classical physics to be a good description). 

Finally, let us point out that the analysis performed in 
this work, and the above quantum CCC, should hold for 
sources with mass to within a few orders of magnitude 
of the Planck mass. Primordial BHs formed in the early 
universe by large density fluctuations could have masses 
in this range. Moreover, it is also plausible that over¬ 
charged configurations with such small masses emerge 
from the gravitational collapse of astrophysical objects, 
acting as seeds for much larger BHs. Our results should 
then apply straightforwardly to these two cases. 


III. CONCLUSIONS 

From the above analysis we can learn two important 
things. First, quantum mechanical effects are perhaps 
able to continuously take us into the classically forbidden 
region of a > 1. This means that even an overcharged 
object, with a charge-to-mass ratio greater than unity, 
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